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Abstract. Mkn 421 was observed during a high flux state for nine nights in April and May 2004 with the fully operational 
High Energy Stereoscopic System (H.E.S.S.) in Namibia. The observations were carried out at zenith angles of 60° -65°, which 
result in an average energy threshold of 1.5 TeV and a collection area reaching 2 km 2 at 10 TeV Roughly 7000 photons from 
Mkn 421 were accumulated with an average gamma-ray rate of 8 photons/min. The overall significance of the detection exceeds 
100 standard deviations. The light-curve of integrated fluxes above 2 TeV shows changes of the diurnal flux up to a factor of 
4.3. For nights of high flux, intra-night variability is detected with a decay time of less than 1 hour. The time averaged energy 
spectrum is curved and is well described by a power-law with a photon index T = 2.1 ± 0.1 sta t ± 0.3 sys and an exponential cutoff 
at E c = 3.1(+0.5 - 0.4) stat ± 0.9 sys TeV and an average integral flux above 2 TeV of 3 Crab flux units. Significant variations of 
the spectral shape are detected with a spectral hardening as the flux increases. Contemporaneous multi-wavelength observations 
at lower energies (X-rays and gamma-rays above » 300 GeV) indicate smaller relative variability amplitudes than seen above 
2 TeV during high flux state observed in April 2004. 
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1. Introduction 

Mkn 421 is a 'BL Lac' type active galactic nucleus. The broad- 
band spectral energy distribution is dominated by non-thermal 
emission that is believed to be produced in a relativistic jet 
pointing towards the observer. The high energy emission of this 
object has been studied by previous observations carried out 
by northern hemisphere ground based Cherenkov telescopes 
(Punch et al. 1992, Aharonian et al. 2002, Krennrich et al. 
2002). Observations of Mkn 421 from the southern hemisphere 
at large zenith angles benefit from considerable increase of the 
collection area at higher energies, which results in a better tem- 
poral resolution at high energies and a better sampling of the 
high energy part of the energy spectrum (Okumura et al. 2002). 

Besides the general interest of understanding the physics of 
the highly relativistic plasma and its interaction with the am- 
bient medium, the proximity of Mkn 421 (z=0.031) makes it 
an interesting target to observe the effect of pair-production of 
gamma-rays with soft (thermal) background photons as part 
of the extragalactic background light. As the collection area 
for large zenith angle observations exceeds two square kilome- 
ters at energies beyond 10 TeV, the observable energy spectrum 
could eventually be extended beyond 20 TeV, which is impor- 
tant to probe the mid-to-far infrared part of the extragalactic 
background light. 

2. Observations and Data analyses 

H.E.S.S. (Hofmann et al. 2003) is an imaging atmospheric 
Cherenkov detector dedicated to the ground based observation 
of gamma-rays at energies above 100 GeV. Situated in Namibia 
(23°16'S 16°30'E), the full array of four telescopes is opera- 
tional since December 2003. Each telescope has a mirror area 
of 107 m 2 and is equipped with a camera consisting of 960 
photomultiplier tubes (Vincent et al. 2003). The system has a 
field of view of 5° and allows to reconstruct the direction of 
individual showers with a precision of better than 0.1°. 

The H.E.S.S. observations reported here were carried out 
for typically 1-2 hours per night from MJD 53107.8 to MJD 
53114.9 (April 12-19, 2004) triggered by an increased level 
of X-ray emission detected by the All-Sky-Monitor (ASM) on- 
board the RXTE satellite and increased activity detected by the 
Whipple Cherenkov telescope (Krawczynski private communi- 
cation). Online analysis of the H.E.S.S. data revealed that the 
source was also active at TeV energies motivating an extension 
of the observational campaign for roughly one week. In the be- 
ginning of May (MJD 53134.8 corresponding to May 8, 2004), 
the H.E.S.S. array participated in a multi-wavelength campaign 
with overlap with pointed X-ray observations with the RXTE 
satellite. 
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The observation mode was employing all four telescopes, 
pointing with 0.5° offset in declination with sign alternating 
from run to run (28 minutes each run). The runs were selected 
according to general quality checks such as absolute value of 
the trigger rate, relative changes in the rate and performance 
of the cameras (number of operational pixels, homogeneity of 
the acceptance). In total, 14.7 hours of good quality data were 
selected for the analyses. The zenith angle of the observations 
ranges from 60.3° at culmination to 65.4° with an average en- 
ergy threshold at 1 .5 TeV. The average size of images and am- 
plitude detected in individual pixels at the given zenith angles 
is well below the saturation limit of the H.E.S.S. cameras and 
read out electronics. 

The data were calibrated following the general procedure 
(Aharonian et al. 2004a) and standard data reduction (image 
cleaning, event reconstruction) was applied (Aharonian et al. 
2005). The cuts used for background rejection were identical to 
the ones used in previous analyses (e.g. Aharonian et al. 2004b, 
2005) with the exception of a loose angular cut 2 < 0.05(°) 2 
(compared to 6 2 < 0.02(°) 2 at smaller zenith angles), where 
8 2 is the squared angular distance between source and recon- 
structed shower direction. The relaxed angular cut compen- 
sates for the degrading angular resolution at low elevations. Its 
value was not optimised to result in a maximum significance 
but rather chosen to ensure a high (> 80 %) and constant accep- 
tance for gamma-rays. The relative energy resolution at the ob- 
served zenith angle range expected from simulation is roughly 
25%. 

For estimating the background, five off regions identical to 
the on region in size and separation from the camera center 
were chosen (Aharonian et al. 2001). The overall signal de- 
tected from Mkn 421 amounts to 7V on = 8978 and </V off > = 
1 357.6 averaged over the five background regions with a signif- 
icance of S = 1 14 cr using the likelihood ratio method derived 
in Li & Ma (1983). 

Finally, in order to reconstruct the energy spectra, collec- 
tion areas were calculated from simulated air showers taking 
the range of observed zenith angles and the offset of the source 
to the camera centre into account. The collection area derived 
for individual events (using linear interpolation between log- 
arithmic energy bins and linear interpolation between zenith 
angles) were used as weights to calculate the flux in bins of 
energy by summing over the events in the on and off regions. 
The effect of spill-over between adjacent energy bins due to the 
instrument's energy resolution was compensated by using col- 
lection areas as a function of reconstructed energy (Mohanty et 
al. 1998, Aharonian et al. 2002). 

Estimates of systematic uncertainties were derived by 
changing the atmospheric transparency in the simulations and 
re-applying the analyses to the data. The small variations ob- 
served in the cosmic ray induced event rate during the observa- 
tions (see also next section) are consistent with smaller varia- 
tions of the atmospheric transparency than assumed in the eval- 
uation of the systematic uncertainties. 
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Fig. 1. Differential energy spectrum of Mkn 421: The curvature is de- 
scribed by a power-law with a photon index of T = 2.1 ± 0.1 stat ± 
0.3 sys with an exponential cutoff at E c = 3.1(+0.5 - 0.4) stat ± 
0.9 sys TeV (solid line) or alternatively with a super-exponential cut- 
off (exp(-(£/£ c ) L6±03 ) at E c = (6.25 ± 0.4 stat ± 0.9 sys ) TeV (dashed 
line). Between 30 and 50 TeV an upper limit with a 99 % confidence 
level is displayed. 

3. Energy spectrum from the Crab nebula at large 
zenith angles 

In order to check the procedure for the reconstruction of energy 
spectra at large zenith angles, a dedicated data-set taken on the 
Crab nebula for 1.8 hours at zenith angles ranging from 57° to 
67° was analysed. 

The fit of a power-law (dN/dE = N Q (E/TeV)- r ) to the re- 
sulting energy spectrum of the Crab nebula yields No = (2.9 + 
0.4 s tat±0.7 sys )x lfr 11 ph/(cm 2 sTeV), F = 2.6±0.1 stat ±0.3 sys , 
X 1 = 7(8 d.o.f.). The result is in good agreement with the re- 
cently published values given by the HEGRA collaboration 
(Aharonian et al. 2004c). This gives confidence in the reliabil- 
ity of the instrument's response and reconstruction technique at 
large zenith angles. 

4. Energy spectra and light-curve from Mkn 421 

The differential energy spectrum for Mkn 421 was calculated in 
the same way as described above for the Crab nebula. The en- 
tire data-set was used to derive a differential energy spectrum. 

The final result for the energy spectrum of Mkn 421 is 
shown in Fig. 1. As the energy spectrum shows clear evidence 
for curvature, the collection areas were calculated iteratively 
until the parameters of the model fit function converged. The 
differential energy spectrum is well sampled between 1 and 
15 TeV. Between 15 and 30 TeV, a signal (corrected for an 
expected spill-over of 3 events) of 2.2cr with N on = 31 and 
(N ff) = 96/5 = 19.2 events is seen. The collection area for this 
energy bin exceeds 2 km 2 after cuts. Between 30 and 50 TeV an 
upper limit with a confidence level of 99 % using the method 
of Feldman & Cousins (1998) was calculated. All energy bins 
are centred on the expected average energy for a given bin as- 
suming a spectrum following a power-law with an exponential 
cutoff. 



The result of a fit to the spectrum assuming a power-law 
with an exponential cutoff (see solid line in Fig. 1), results 
in N = (1.55 ± 0.08 stat ± 0.4 sys ) x 10~ 10 ph/(cm 2 s TeV), 
F = 2.1±0.1 stat ±0.3 sys , E c = 3.1(+0.5 -0.4) stat +0.9 sys TeV, and 
X 1 = 16(10 d.o.f.). Other model fits using more complicated 
functions which have been commonly used in the past to de- 
scribe a curved spectrum were tested: a parabola in \og(dN/ dE) 
vs. log(£') gives x 2 = 12.6(9 d.o.f.) and a super-exponential 
cutoff^ 2 = 1 1.7(9 d.o.f). The super-exponential cutoff fit with 
dN/dE oc E- T x exp(-(E / E c ) a ) with a = 1.6 ± 0.3 (systematic 
uncertainty negligible) is included in Fig. 1 as a dashed line. 
According to the F-test, the probability of a chance improve- 
ment of the^ 2 for the super-exponential cutoff is 6.4 % which 
is not sufficiently low to give a preference to this fit function 
(the same holds true for the parabolic fit function). 

The cutoff energy E c = 3.1(+0.5 - 0.4) stat + 0.9 sys TeV 
found in the observations described here is fully consistent with 
previous observations in the years 2000 and 2001 by HEGRA 
of 3.6(+0.4 -0.3) sta t(+0.9 -0.8) sys TeV (Aharonian et al. 2002) 
and VERITAS (4.3 ± 0.3) TeV (Krennrich et al. 2002). 

The diurnal integral fluxes above 2 TeV were calculated by 
taking the sum over the inverse collection area for the events 
with energies above the chosen value. With the threshold en- 
ergy raised to 2 TeV, this gives an estimate of the flux which 
is independent of the variation of the actual energy threshold 
within a night, at the expense of loosing count statistics. As 
an additional benefit, the dependence on changes of the spec- 
tral shape are negligible as the collection area is not strongly 
changing with energy for energies well above the threshold. 
The resulting light-curve is shown in Fig. 2. 

The light-curve exhibits night-by-night variability in the 
first observing week and resumes a lower flux in the later obser- 
vation. The peak diurnal average flux reaches a value of about 
5 times the flux observed from the Crab nebula. The corre- 
sponding gamma-ray rate is sufficiently high to probe intra- 
night variability of the high energy end of the spectrum with 
unprecedented accuracy. 

During the observations of MJD 53113.8-53113.9 (April 
18) and 53114.8-53114.9 (April 19), significant variations of 
the flux within these nights are detected. The hypothesis of a 
constant flux during these nights results in x 2 = 40.6(7 d.o.f.) 
and x 2 = 28(7 d.o.f.) respectively. The intra-night variations 
as seen during MJD 531 13.8-531 13.9 (April 18) using bins of 
14 min width are shown as an inlay in Fig. 2. In order to exclude 
variations of the detector's response or changes in the atmo- 
sphere to be responsible for the observed variations, the post- 
cut (after applying the image shape cuts) cosmic ray rate has 
been checked for variability which is smaller than 2 % in rel- 
ative root mean square (RMS) during this night. The observed 
significant intra-night variability in the light-curve suggests a 
decay time of less than 1 hour. 

The light-curve is not corrected for possible variations of 
the atmospheric transparency. The atmospheric transmissivity 
between the position of the air shower and the detector can 
be probed by the rate of detected cosmic ray events, as the 
Cherenkov light traverses the same aerosol layer which is be- 
lieved to dominate possible temporal variations of the atmo- 
sphere's transparency. During the observations, the average 
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Fig. 2. The light-curve for Mkn 421 in April and May with diurnal 
integral flux values above 2 TeV. For comparison, the horizontal line 
indicates the flux observed from the Crab nebula. The inlay shows the 
light-curve for the night of MJD 53113 (marked with a circle and an 
arrow) with the strongest indication for intra-night variability in 14 
min bins (0.01 d). 

cosmic ray rates for individual runs varied between 113 and 
143 Hz with an average of 132 Hz and RMS of 13 Hz. The 
post-cut rate shows even less variation with 8 % relative RMS. 
The observed variability of Mkn 421 is therefore clearly asso- 
ciated with the source and not a consequence of temporal vari- 
ations of the detector response or atmospheric transparency. In 
principle, a correction of the measured gamma-ray flux could 
be derived and applied to the light-curve. This was not done 
as the effect is small (< 10%) in comparison to the observed 
variability. 

In order to study variations of the spectral shape in vari- 
ous flux states, diurnal spectra were calculated. Given the cor- 
relation between the cutoff energy E c and the photon index T 
derived from the fit of a power-law with an exponential cutoff 
(correlation coefficient -0.955 between E^ 1 and photon index 
T), the power-law index was kept constant (F = 2) for the fits 
applied to spectra obtained for individual nights. The fit was 
applied to a fixed energy range from 1 to 10 TeV. As a result, 
a hardening of the spectrum or an increase of the cutoff energy 
is seen clearly in the correlation of the cutoff energy E c and 
the integral flux above 2 TeV in Fig. 3. A similar result was 
obtained when keeping the cutoff energy fixed and letting the 
photon index vary freely. 

5. Multiwavelength observations in April 2004 

During April 2004, a coordinated multi-wavelength campaign 
monitored the activity of Mkn 421 in radio, optical, X-ray, and 
gamma-rays (Cui et al. 2004). The source was seen to be active 
in X-rays where observations with the array of proportional 
chamber units (PCU) onboard the RXTE satellite were per- 
formed. These observations were not simultaneous with the ob- 
servations with the H.E.S.S. array, but by combining the PCU 
with the ASM data a good temporal coverage overlapping with 
the H.E.S.S. observations can be achieved. An average of 4 
counts/s were detected by the ASM during the first weeks of 
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Fig. 3. The cutoff energy as determined from a fit with a fixed power- 
law photon index T = 2 versus integrated flux measured for E > 
2 TeV. Note, that the point with the lowest flux is from the May obser- 
vations. 

April. This is sufficiently high to probe the activity of Mkn 421 
during individual ASM pointings (dwells). 

The Whipple 10 m Cherenkov telescope was observing 
Mkn 421 simultaneously with the RXTE pointings which were 
generally starting within a few hours after the H.E.S.S. obser- 
vations took place. First preliminary results of this campaign 
have been presented by Cui et al. (2004). 

In Fig. 4, the different observations are combined such that 
the observed flux (or count-rate) is normalized to the average 
flux (or count-rate) during the time between MJD 53107 and 
53116. The preliminary Whipple light-curve is derived from 
the count-rate which is not corrected for different zenith angles 
of observations. Generally, the Whipple observations were car- 
ried out at small zenith angles (Cui et al. 2004) and according to 
previous observations, the energy threshold is estimated to be 
« 400 GeV (Krennrich et al. 2002). The ASM light-curve and 
the 1 cr uncertainty band is obtained by calculating the sliding 
average over five dwells. 

Even though a detailed correlation analysis is difficult be- 
cause of a large fraction of the data not being simultaneous, it 
is interesting to note that the strongest variability occurs at the 
highest energies, with a maximum relative amplitude of the ob- 
served flux of fmax/^min = 4.3 + 0.5 as compared to 1.7 + 0.2 
derived from the Whipple count rates and 2.05 + 0.02 from 
the RXTE PCU observations. As a measure of the variability 
the relative RMS of the measurement indicates a variability of 
(17 + 6) % seen by the Whipple instrument and (26 ± 1) % 
for the X-ray data taken with the PCU detectors, whereas for 
H.E.S.S. the relative RMS value is (51 + 8) %. 

6. Conclusion 

The results on the time averaged energy spectrum of Mkn 42 1 
presented here confirm previous observations of the existence 
of a cutoff at E c = 3.1(+0.5 - 0.4) stat + 0.9 sys TeV which is 
lower than the cutoff energy of 6.2 ± 0.4 stat (+2.9 - 1.5) sys TeV 
observed from Mkn 501 in 1997 (Aharonian et al. 1999). Given 
the similar red shift of these two objects, this would imply that 
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Fig. 4. The relative changes of the flux with respect to the average ob- 
served during the multi-wavelength campaign in April 2004. Upper 
panel: H.E.S.S. F(E>2 TeV) (filled circles), Whipple (crosses, prelim- 
inary count rate without correction for varying zenith angles). Lower 
panel: PCU (open triangles, count-rate), ASM (1 cr error band, see 
text for details). Whipple and RXTE PCU data taken from Cui et al. 
(2004). 



the cutoff in the Mkn 421 spectrum is intrinsic to the source 
and not due to absorption. With the observations carried out at 
large zenith angles and increased collection area, a signal at the 
level of 2.2 cr (corrected for the effect of spill-over) was de- 
tected between 15 and 30 TeV. This detection in conjunction 
with the observation of a cutoff in the energy spectrum intrin- 
sic to the source is important for constraining the extragalactic 
background light at wavelengths beyond 10 can. Furthermore, 
the energy spectrum was observed to become harder as the in- 
tegral flux increases. It is not possible to discern whether the 
hardening is a consequence of an increase of the cutoff energy 
from 1 .5 to 3.5 TeV or a change in the power-law index as these 
two parameters are highly correlated. 

The relative amplitude and variance of the variability ob- 
served above 2 TeV is significantly higher than observed at 
lower energies during the same activity period by the Whipple 
telescope and the RXTE pointed X-ray detector. The sparse but 
similar sampling of the light curve among the different obser- 
vations (H.E.S.S., Whipple, and RXTE PCU) does not intro- 
duce a bias which could explain the observed difference in the 
variability. 

The larger variability seen at multi-TeV energies compared 
to energies below TeV is consistent with the observed spectral 
changes, where the relative increase of the flux at higher en- 
ergies causes a hardening of the observed spectra. Assuming 
the correlation between integral flux and cut-off energy indi- 
cated in Fig. 3 and using the integral flux above 2 TeV as mea- 
sured with H.E.S.S., the expected light curve above 400 GeV 



was calculated. As expected, the relative RMS scatter of the 
lower energy light curve is roughly half of the value observed 
above 2 TeV well consistent with the observed variability in the 
Whipple light curve. 
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